Spin Gap in a Doped Kondo Chain 
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H = -tY^ c\„Cj„ + h.C. + JK^S^ ■ Si + Jn^S^ ■ Sj. 

{ij)a i {ij) 
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Here c^ is a creation operator for an electron of spin a 



J 1, 



is the conduction electron 
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We show that the Kondo chain away from half-filling has a spin gap upon the introduction 
of an additional direct Heisenberg coupling between localized spins. This is understood in the 
weak-Kondo-coupling limit of the Heisenberg-Kondo lattice model by bosonization and in the strong- 
coupling limit by a mapping to a modified t-J model. Only for certain ranges of filling and Heisenberg 
coupling does the spin gap phase extend from weak to strong coupling. 

PACS numbers: 75.20.Hr, 75. 30. Mb, 74.20.-z 

The Kondo lattice model (KLM) [Eq. (0) with Jh = 0] 
is often considered as a model for heavy fermion materi- 
als. The one dimensional (ID) version has recently been 
argued to be very relevant to a theoretical approach to 
high-Tc superconductivity based on fluctuating hole-rich 
stripes [Q. It has been rather clearly established that, 
at half-filling, the ID KLM has a gap to both spin and 
charge excitations m. On the other hand, its behavior 
away from half-filling remains controversial. Of particu- 
lar interest is the possibility of a phase with a spin gap 
but no charge gap. In this case, the charge density wave 
and superconducting susceptibilities are significantly en- 
hanced, and for effectively attractive interactions the sys- 
tem is expected to become superconducting upon the in- 
troduction of any inter-chain hopping. 

A large Kondo coupling forces all conduction electrons 
to form on-site singlets with localized spins. The un- 
paired localized spins then effectively hop on the back- 
ground of singlets via a reduced hopping term t/2. Ad- 
ditional t^ / Jk interactions then lead to incomplete fer- 
romagnetism B . At weaker Kondo coupling a Luttinger 
liquid (LL) phase has been proposed and supported, to 
some extent, by density matrix renormalization group 
(DMRG) work |,|. 

Part of the difficulty in understanding this model is 
that it cannot be bosonized in a very direct manner. Two 
approaches have been proposed. In one of them [|j a spe- 
cial, "Toulouse" , limit with a highly anisotropic Kondo 
interaction is bosonized, resulting in a spatially modu- 
lated sine-Gordon model and leading to the prediction 
of a spin gap (SG) phase with a spin gap proportional 
to the g power of the perpendicular part of the Kondo 
coupling. 

An alternative approach can be applied to a more gen- 
eral Heisenberg-Kondo lattice model (HKLM) which has 
a direct Heisenberg coupling Jh between the localized 
spins as well as the usual Kondo interaction: 



at site i, Si = ^^^ c:„tcr„0C,^ 

spin operator. Si is the localized spin-i operator, and {ij) 
denotes nearest neighbors. This model can be bosonized 
MM in the limit Jk ^ t, Jh in which only the low en- 
ergy degrees of freedom of the Heisenberg chain and the 
conducting free electron chain are kept. This approach 
predicts JS) a SG phase. Another interesting limit of the 
generafized model H,^ is Jk ^ t and Jk ^ Jh ^ t^ / Jk 
in which it becomes equivalent to a t-J model with an 
effective electron density of ncff — \1 ~ n\, where n is the 
conduction electron density. Building on the rather well 
understood phase diagram of the t-J model we learn that 
LL and phase-separated (PS) regimes exist as well as a 
SG phase for some range of n near half-filling. Thus we 
see that there must be at least four phases in the HKLM. 

We note that the usual KLM contains an effective 
RKKY interaction between the localized spins, but this 
would normally be of longer range and weaker than our 
Jh. Such a short-range exchange interaction is believed 
to be present, although very small, in certain quasi-one- 
dimensional organic conductors containing Cu^+ ions, 
Cu(tatbp)I and Cu(pc)I Q. 

The purpose of the present work is to check the an- 
alytical predictions based on bosonization and the large 
Jk limit using the DMRG numerical method. In partic- 
ular, we establish the existence of the SG phase over a 
significant region of parameters. We first briefly review 
the large Jk limit and the bosonization approach of Ref. 
Q, clarifying some points about the large Jk Hmit. 

When Jk » t and Jk > Jh » t^ / Jk, the HKLM 
reduces to the t-J model P,B| except that the Coulomb 
term niUj usually included in the t-J model is absent: 
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with the rc-dcfinitions t — *■ t/2 and n — > |1 — n|, and Jh 
playing the role of J. The no-double-occupancy restric- 
tion in the t-J model is automatically met by the HKLM 
in the limit of large Jk since the effectively- hopping spins 
have the same restriction. Unlike in the t-J model, the 
Heisenberg coupling exists not only between the unpaired 
spins; however, the connection with the t-J model is seen 
to be complete upon noting that the Heisenberg cou- 
pling is ineffectual for localized spins which are already 
strongly coupled with the conduction electrons in an on- 
site singlet pair. 

Numerical work on the one-dimensional t-J model |0] 
show a spin gap phase for low electron density and 
2 < J/t < 3 (Fig. |l|). This SG phase adjusts somewhat 
upon dropping its Coulomb term [^, as can be seen by 
considering the case of two electrons in the t-J model 
following Ref. [O. Using the ansatz ground state 
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where a is a constant to be determined and b^^ = 

'^ii^il'^i+n i ~'^ii^i+n t) ^^ ^ singlet statc creation opera- 
tor giving two electrons n lattice sites apart, we find that 
the boundary between the LL and SG phases in the limit 
of zero filling in the t-J model is shifted upon dropping 
its Coulomb term from J — 2t to J = 8t/3, and a rough 
overestimate [see Fig. |l| of that between the SG and PS 
phases is shifted from J = 3.2t to J = 6.3i 0. Using this 
as a guide, we can choose values of Jh and filling which 
allow us to find a spin gap in the strong-Kondo-coupling 
limit of the HKLM. 

The phase diagram of the HKLM is three dimensional, 
with parameters Jk/^, Ja/t, and the filling n/2. Since 
filling is the most difficult parameter to vary in our nu- 
merical DMRG study, we choose a particular value n = ^ 
to ensure that in the large Jk limit we can pass through 
the SG region known from the t-J model by varying Jh. 
At large Jk we expect, as we increase Jh from zero, to 
find first a small ferromagnetic (FM) region until the 
Heisenberg term becomes comparable to the 0{t? / J^i) 
fluctuations responsible for ferromagnetism, then a LL 
region, followed by a SG region, and finally a PS region. 

Thinking of the HKLM as a spin-i chain coupled to a 
conduction electron chain by a weak interchain (Kondo) 
coupling, we proceed via non-Abelian bosonization, writ- 
ing 



(4) 



Sj ~ pI + Pr + ai{-iytY{agi) 
for the spin chain, and 

a2[e^"'^='tr(a92)e'^^'^=-Fh.. 



for the conduction electrons. {ai,2 are constants, and 
we employ the convention of using the variable x for 
the continuum representation of lattice position j.) The 
SU(2) currents p^'^{x) have scaling dimension 1, while 
the SU(2) matrix fields 171,2(2^) and the charge boson op- 
erator e'v27r0c(a:) jjg^yg dimension i. At half- filling the rel- 
evant (scaling dimension |) product of alternating terms 
produces a gap, Ag oc J^. Away from half- filling, below 
the length scale l/(7r — 2fci?), only the marginal interac- 
tion 
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survives, similar to the case of the spin- 2 spin ladder with 
zig-zag rungs |q| . Under renormalization group flow, this 
interaction flows to zero coupling for FM Jk < 0, giv- 
ing no gap, and to strong coupling for antiferromagnetic 
Jk > 0, giving a spin gap 
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where c is a positive constant, v = vi -\- V2 0, and the 
two chains' different spin-wave velocities are vi = tt Jh/2 
and V2 = 2t sin(7rn/2). 

AU of the DMRG calculations for the HKLM reported 
here were done using the finite system method p3| . Spin 
gaps are obtained from reliable extrapolations to infinite 
chain length from several different lengths commensu- 
rate with the filling n — -^ admitting total spin pro- 
jection S^ = 0, namely L = 16, 32, 48, and 64. The 
calculations were done on an IBM RS/6000 workstation 
equipped with 64MB RAM; typical processing times are 
on the order of 20 hours of CPU time per ground state 
calculation for a chain of length 64 keeping 180 states 
in two finite system method sweeps, although such ac- 
curacy was not required for all points calculated. As a 
general rule, the discarded weights (sums of the density 
matrix eigenvalues neglected in the basis truncation at 
each step) in this case are usually on the order of 10~® 
for large Kondo coupling Jk and up to 10~^ for small Jk 
or near phase boundaries. 

Our DMRG calculations verify that there are points 
in the phase diagram which have a spin gap, and points 
which do not. Where these are close enough together, 
we can estimate the location of the phase boundary. 
Based on our results, we find this boundary along the 
line Jk = lOt lies between Jh = 1.25t and 1.65i, and 
along the line Jh = 1.25i it lies between Jk ~ 5t and 
Jk = lot (close to Jk — 7.5i). As shown in Fig. ||, the 
phase boundary along the line Jk = 2t is near Jh = 0.6i. 

A test of the bosonization prediction of Eq. (0) is de- 
picted in Fig. pi For large Jk the exponential prediction 
is not applicable and in fact the gap goes through a max- 
imum in the approach to the large Jk limit which maps 
onto a modified t-J model. The resulting approximate 
linearity in the range t < Jk < 3t is good verification of 



the bosonization prediction. The slope is to be compared 
with —cv\ the ratio of the slopes for the two values of Jh 
is predicted to be 0.86, while it is numerically determined 
to be 0.83 ± 10%. The value of c determined from the 
DMRG calculations is 1.1 ± 10%. 

We are thus led to conjecture a phase diagram for the 
HKLM, near half-filling, like the one shown in Fig. 0. 
We can regard the large Jk region of this diagram as 
being rather firmly established by prior work on the t- 
J model. We have also assumed the correctness of the 
bosonization prediction of a spin gap at small Jk (for 
all non-zero Jh), encouraged by the above DMRG ver- 
ification. This DMRG work also suggests that the SG 
phase at large Jk is connected to the one at small Jk. 
The boundary of the FM phase at a finite J^ has also 
been established by previous work. We expect this FM 
phase to persist, for large enough Jk, up to some finite 
Jh where the predominant interaction between unpaired 
localized spins becomes antiferromagnetic. Possibly the 
most interesting phase boundary is the one between LL 
and SG phases. This could extend all the way to Jk = 
as drawn in Fig. 0, in which case the SG phase would not 
occur in the pure Kondo lattice model. Alternatively it 
could intercept the Jk axis at a finite value, J^ < J^. 
Finally, it might intercept the boundary of the FM region 
at finite Jh in which case the pure KLM would make a 
transition directly from the FM phase into the SG phase 
without passing through the LL phase. In the latter two 
scenarios, the spin gap in the pure KLM should turn on 
below Jj^ or J^ , pass through a maximum, and then van- 
ish at Jk — > where the model becomes non-interacting. 
If this maximum spin gap is not too small (greater than 
a few percent of t) it should be observable by DMRG. 
An important parameter for the spin excitations in ei- 
ther LL or SG phase is the spin- wave velocity, Vs- In the 
pure KLM, we expect this to go to zero at the bound- 
ary of the FM phase, pass through a maximum as Jk 
decreases, and then vanish as Jk — ^ 0. In the second 
scenario, it should remain finite at J^, where the spin 
gap turns on. We expect that the spin gap will not ex- 
ceed Vs in order of magnitude. In regions where Vs is 
very small (which may include the entire Jk axis below 
Jj^) the LL (or SG) behavior only appears at very low 
energies (and long lengths) and hence is difhcult to study 
numerically. A non-zero Jh enhances Ug, making it eas- 
ier to determine the phase diagram. So far, we have not 
found a spin gap in the pure KLM, although we have not 
exhaustively studied all values of Jk and n. For n — ^, 
we see evidence of a FM/LL phase boundary between 
Jk/^ — 1 and 2; no spin gap is seen down to Jk/^ — j- 
This is consistent with what was found earlier by DMRG 
in a somewhat different parameter regime H . For lower 
densities, the FM/LL boundary moves closer to Jk = 0, 
since in the limit of n -^ the ground state of the KLM 
has been rigorously shown to be FM uM, closing the 
window of opportunity for a SG phase in the KLM. At 



densities lower than about n = |, we know from work on 
the t-J model that the SG phase doesn't exist at large 
Jk. Nevertheless, assuming the veracity of the bosoniza- 
tion approach, we might expect the SG phase to exist at 
non-zero Jh and smaller Jk. 

The above bosonization approach implies that the SG 
phase is a charge-only LL. This conclusion is supported 
by the fact that we find numerically that the charge- 
transfer gap Act{n,L) = EL{nL + 2) + EL{nL — 2) — 
2EL{nL), where El{Q) is the ground state energy of the 
length-L chain in the sector of charge Q and S^ = 0, is 
zero in the thermodynamic limit in this phase. 

A conformal field theoretic finite-size scaling analy- 
sis predicts that in a charge-only LL, the coefficient of 
Q^ in the energy of the ground state in the sector of 
charge Q should be ^ • j-Kvi -f • ■ ■ where • ■ • repre- 



sents terms of higher order in 1/L, including 1/LlnL 
terms as well, and Vc is the charge velocity. Selecting 
a point (Jk, Jh) = (3i, 1.65i) at which the spin gap is 
reasonably large we determined El{Q) by the DMRG 
method for wide ranges of Q for i = 32 and 64, and 
find that the finite-size scaling is well verified; both chain 
lengths studied give Vc/{v^'iTrR'^) = 0.082 ± 3% where 
Vp ^ V2 ~ 2tsm{TTn/2) is the charge, or Fermi, velocity 
in the non-interacting case. 

In conclusion, we have verified that the SG phase of 
the HKLM exists over an extended region of parameter 
space interpolating between the large Jk region where it 
follows from previous work on the t-J model to the weak 
Jk region where it is predicted by bosonization. We have 
verified the expected exponential dependence of the gap 
on v/ Jk in the weak coupling limit. So far, we have found 
no evidence for a spin gap in the pure KLM but further 
numerical work is required (and is in progress) to settle 
this question one way or the other. 
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FIG. 3. Testing the small- Jk bosonization prediction for 
two values of Jh by plotting ln(As/i) versus t/ .Jk- 
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FIG. 1. A numerically-computed phase diagram of the 

one- dimensional t-J model (taken from Ref. |[ll|). 




FIG. 2. The spin gap in the HKLM along Jk = 2i as a 
function of Jh shows a transition near Jh = 0.6i. Error bars 
are approximate measures of how much we could reasonably 
expect the values to change by keeping all states. 
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FIG. 4. A possible phase diagram of the one-dimensional 
HKLM at n = |. Filled circles are points at which a finite 
spin gap was measured by DMRG, open circles are gapless 
points (or have spin gaps less than 0.5% of t), and shaded cir- 
cles are too close to a phase boundary to determine the gap. 
For fillings further from n = 1, the spin gap phase at large Jk 
shrinks and eventually closes off, while remaining at smaller 
Jk. 



